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'qq" ■ ABSTRACT 

We present a simple set of kinematic criteria that can distinguish between galaxies dominated by 
ordered rotational motion and those involved in major merger events. Our criteria are based on the 
CM ' dynamics of the warm ionized gas (as traced by Ha) within galaxies, making this analysis accessible 

to high-redshift systems, whose kinematics are primarily traceable through emission features. Using 
the method of kinemetry (developed by Krajnovic and coworkers), we quantify asymmetries in both 
pLn \ the velocity and velocity dispersion maps of the warm gas, and the resulting criteria enable us to 

empirically differentiate between non- merging and merging systems at high redshift. We apply these 
criteria to 11 of our best-studied rest-frame UV/optical-selected z ~ 2 galaxies for which we have near- 
infrared integral field spectroscopic data from SINFONI on the VLT. Of these 11 systems, we find 
■ that > 50% have kinematics consistent with a single rotating disk interpretation, while the remaining 

systems are more likely undergoing major mergers. This result, combined with the short formation 
timescales of these systems, provides evidence that rapid, smooth accretion of gas plays a significant 
. role in galaxy formation at high redshift. 

Subject headings: methods: data analysis — techniques: spectroscopic — galaxies: evolution — galax- 
ies: high-redshift — galaxies: interactions — galaxies: kinematics and dynamics 



1. INTRODUCTION 

In recent years, deep observations of galaxies at 
z ~ 2 — 3 have provided detailed insight into the growth 
of structure and galaxy evolution in the early Universe. 
High-resolution broad-band imaging and long-slit spec- 
troscopic surveys have shown the population at this red- 
shift to be rapidly evolving and diverse. This is a reflec- 
tion of the significant growth of galaxies that occurs at 
this epoch; at z ~ 2, both the cosmic star formation rate 
and the luminous quasar space density are at their peak 
(e.g. Fan et al. 2001; Chapman et al. 2005; Hopkins & 
Beacom 2006). Correspondingly, the stellar mass den- 
sity in galaxies increases from ~ 15% its current value at 
z ~ 3 to 50 — 75% its current value at z ~ 1 (e.g. Dickin- 
son et al. 2003; Fontana et al. 2003; Rudnick et al. 2003, 
2006), making this the era when much of the assembly 
of massive galaxies occurs. Systems at this redshift con- 
sequently show a wide range in properties, with nuclear 
activities varying from negligible to active to powerful 
QSOs, star formation rates varying from less than 1 M 
yr~ 4 (e.g. passively evolving BzK-selected objects and 
quiescent DRGs; Cimatti et al. 2004; Daddi et al. 2005; 
Labbe et al. 2005; Kriek et al. 2006b; Wuyts et al. 2007) 
to over 10 3 M yr^ 1 in sub-mm-selected galaxies (SMGs; 
e.g. Blain et al. 2002; Smail et al. 2002; Tacconi et al. 
2006), and correspondingly large variations in morphol- 
ogy, stellar populations, excitation properties, and dust 
content (e.g. Reddy et al. 2005; Papovich et al. 2006; 
Kriek et al. 2006a). 

Furthermore, the resolving power of 8-10m class tele- 
scopes and of millimeter interferometric arrays reveals 
velocity gradients within many of these systems (Erb 
et al. 2003; Genzel et al. 2003). By coupling such tele- 
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scopes with high-resolution integral field spectrographs, 
we are now able, for the first time, to resolve the dy- 
namic structures and internal processes at work within 
massive galaxies during their critical stages of evolution 
(Forster Schreiber et al. 2006; Genzel et al. 2006; Wright 
et al. 2007; Law et al. 2007, see also e.g. Flores et al. 
2004; Puech et al. 2006; Swinbank et al. 2006). With sur- 
veys of the spatially- resolved kinematics of various high-z 
populations (Bouche et al. 2007), we can now begin to 
understand the forces driving such rapid and intense evo- 
lution as well as the role of secular evolution and major 
mergers in these processes. 

Differentiating between systems in ordered rotation 
and those undergoing major merger events has signifi- 
cant ramifications in understanding the evolution of both 
the baryons and the underlying dark matter distribu- 
tions. Kinematic measurements of a system's baryonic 
component, combined with basic assumptions about the 
morphology of the system, enables a detailed probe of 
the mass and angular momentum of dark matter halos 
at z ~ 2 and of the interaction between these halos and 
their baryons (Forster Schreiber et al. 2006; Bouche et al. 
2007). Additionally, the baryons themselves can also 
constrain formation and evolution scenarios, by probing 
whether the active star formation seen in these objects 
is triggered by major mergers (e.g. Hopkins et al. 2006) 
or by smooth accretion (Birnboim et al. 2007). In sev- 
eral well-resolved systems at z ~ 2, the young age of 
the stellar population (~ 500 — 1000 Myr), when com- 
bined with the high star formation rate (up to ~ 200 
M© yr" 1 ), suggests extremely rapid (< 1 Gyr) forma- 
tion (Forster Schreiber et al. 2006; Genzel et al. 2006). 
This is most surprising since the dynamics of these sys- 
tems qualitatively appear to be consistent with no recent 
major merger events, thus indicating a rapid and intense, 
but still smooth, mass accretion mechanism. 

However, a quantitative and definitive understand- 
ing of the structure of these high-z systems is compli- 
cated by the limited spatial resolution attainable at this 
rcdshift and the lower signal-to-noise (S/N) associated 
with these faint objects. Existing prescriptions for mea- 
suring galaxy morphology use surface brightness infor- 
mation and rely solely on characteristics of the broad- 
band emission (CAS: Conselice 2003; Gini/M20: Lotz 
et al. 2004; Sersic fitting: Ravindranath et al. e.g. 2004; 
Cresci et al. e.g. 2006). At lower redshifts, these tech- 
niques have been shown to effectively distinguish dis- 
parate populations over a wide range of resolutions and 
S/N (Lotz et al. 2004). At z > 2, however, the sit- 
uation is more complicated, since optical observations 
probe the rest-frame ultraviolet morphology, which is 
strongly affected by extinction and by the light from 
massive stars in star-forming regions. Near-infrared ob- 
servations, which correspond to rest-frame optical emis- 
sion at z ~ 1—4, are a significant improvement, but 
high-resolution near-infrared imaging is still quite obser- 
vationally expensive. To reliably probe the dynamical 
state of a system, a promising alternative comes from 
integral-field spectroscopic (IFS) observations, yielding 
spatially-resolved kinematic information. While this ap- 
proach is also still observationally expensive, it has the 
advantage of directly probing the system's dynamics and 
total enclosed mass. 

Given the detailed kinematic information available 



with such observations, a technique complementary to 
morphological methods can be developed to fully exploit 
the measured two-dimensional velocity structure. Exist- 
ing IFS observations of Ha emission in a few cases at 
z ~ 2 reveal dynamics suggestive of either the "spider 
diagram" structure found in local disk galaxies or of the 
complex structures found in mergers (Forster Schreiber 
et al. 2006). If the observed internal dynamics of the 
warm gas do in fact reflect those of the system, then it 
appears possible to distinguish mergers and non-mergers 
using current data. Indeed, simulated observations with 
current technology have predicted that the kinematic 
differences between non-merging and merging systems 
should be qualitatively visible in IFS data (Law et al. 
2006). 

In this paper, we present a scheme for discerning be- 
tween merging and non-merging galaxies, based on their 
emission-line kinematic properties and also on the dis- 
tribution of the stellar continuum intensity, which, with 
integral-field observations, can be unambiguously sepa- 
rated from the emission lines. We use kincmetry (Kra- 
jnovic et al. 2006) to quantify asymmetries in the velocity 
and velocity dispersion maps, enabling us to differentiate 
a system in regular, ordered rotation from one disturbed 
by the complex dynamics of a major merger, even at 
the ~ 4 kpc spatial resolution typical of seeing-limited 
(FWHM ~ 0'.'5) observations of z ~ 2 systems. In Sec- 
tion 2, we describe this method, and in Section 3, we il- 
lustrate the power of our technique with a number of tem- 
plate galaxies, drawn from observations of local systems 
and from simulations (from Daigle et al. 2006; Colina 
et al. 2005; Naab et al. 2007). Section 4 applies this clas- 
sification scheme to the well-resolved z ~ 2 systems from 
the Spectroscopic Imaging survey in the Near-infrared 
with SINFONI (SINS; see Forster Schreiber et al. 2006; 
Genzel et al. 2006; Bouche et al. 2007; Cresci et al. in 
prep). Section 5 discusses these results and addresses the 
inherent assumptions in our method, including the use of 
Ha kinematics as a tracer of dynamics at high redshift. 
Finally, we summarize our conclusions in Section 6. 

Throughout this paper, we compare the two classes 
of systems that we aim to distinguish: those with re- 
cent major merger events (mass ratio < 3 : 1) and those 
without. For simplicity, we refer to these two classes 
as "mergers" and "disks" respectively, but the obvi- 
ous caveats of such nomenclature are worth mentioning. 
Our analyses focus mainly on the gaseous component 
in galaxies, which relaxes quickly into a flat, disk-like 
configuration, e.g. in late-stage mergers (e.g. Barnes & 
Hernquist 1996; Mihos & Hernquist 1996; Tacconi et al. 
1999; Naab et al. 2006). The naming scheme adopted 
here does not reflect an inability to distinguish late-stage 
mergers from non-merging systems; in fact, Section 3 il- 
lustrates that the gas motions in even late-stage mergers 
are still sufficiently disturbed that they can be clearly 
identified as such (see also Barnes 2002; Arribas & Col- 
ina 2003). In the context of this paper, these systems 
are never referred to as "disks," despite their probable 
geometry; that term is instead reserved for galaxies that 
show no sign of a recent major merger. 

We assume a A-dominated cosmology with H = 70 km 
s _1 Mpc -1 , O m = 0.3, and Cl\ = 0.7. For this cosmology, 
1" corresponds to w 8.2 kpc at z = 2.2. 
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Fig. 1. — From left to right: Continuum intensity in rest-frame ii-band, intensity of the Ha line emission, and velocity and velocity 
dispersion of the Ha-emitting component, for a simulated ideal disk (top) and for the high-redshift galaxy BzK-6004 observed as part of 
the SINS program (bottom). In the ideal disk, the star formation (Ha intensity) follows the mass distribution (Continuum), whereas in 
the observed z ~ 2 system, a significant off-center star-forming region is seen. Overplotted on the velocity maps are isovelocity contours, 
whose concave curvature on either side of the rotation axis displays the so-called "spider diagram" structure characteristic of rotational 
motion (van dcr Kruit & Allen 1978). Also overplotted on the velocity and dispersion maps are sample ellipses from the expansion with 
kinemetry. Along a kinemetry ellipse in the velocity map, the velocity varies as the cosine of the azimuthal angle. Along an ellipse in the 
velocity dispersion map, the dispersion is approximately constant with angle. 



2. METHOD 

To determine whether a particular observed system is a 
disk or a merger, we use two main criteria: the symmetry 
of the velocity field of the warm gas, and the symmetry 
of the velocity dispersion field of the warm gas. An ideal 
rotating disk in equilibrium is expected to have an or- 
dered velocity field, described by the so-called "spider 
diagram" structure, and a centrally-peaked velocity dis- 
persion field (Figure 1, see also van der Kruit & Allen 
1978). 

Likewise, such a disk is also expected to have the reg- 
ular and centrally-peaked continuum distribution char- 
acteristic of exponential disks. Indeed, this feature is 
the basis for low-redshift morphological classification 
schemes (e.g. Conselice 2003; Lotz et al. 2004; Cresci 
et al. 2006). However, in our IFS observations, the high- 
est S/N is usually obtained from the Ha emission from 
the warm, star-forming gas, rather than from the un- 
derlying stellar continuum whose light is dispersed over 
many more spectral pixels. In Forster Schreiber et al. 
(2006), for example, the emission features have a typical 
S/N of ~ 25, whereas the stellar continuum is often only 
detected in a small part of a given system, with typical 
S/N < 10 in the brightest region. We therefore perform 
only the simplest analysis of the distribution of the stel- 
lar component in these systems, using this information 
to supplement and inform the detailed analysis possible 
with the high-S/N emission line velocity and velocity dis- 
persion data. 

In our analysis, we do not include constraints based 
on the intensity distribution of the emission lines. This 
tracer of the location and power of star-forming regions 



is often clumpy and asymmetric in even the most kine- 
matically regular disks (Figure 1; see also Daigle et al. 
2006 for local examples) and, consequently, reveals lit- 
tle about the mass distribution and dynamical state of 
the system. Rather, it is the kinematics of this gas, and 
not its spatial distribution, that reflects the dynamical 
state of the system and therefore forms the basis for our 
analysis. 

2.1. Quantifying Symmetries with Kinemetry 

The symmetries in kinematic fields can be measured 
via the kinemetry method developed and described in 
detail by Krajnovic et al. (2006). Briefly, kinemetry is 
an extension of surface photometry to the higher-order 
moments of the velocity distribution. The procedure op- 
erates by first describing the data by a series of concen- 
tric ellipses of increasing major axis length, as defined by 
the system center, position angle, and inclination. The 
latter two parameters can cither be determined a pri- 
ori and used as inputs or can be measured functions of 
semi-major axis length as a first step in the kinemetric 
analysis. Along each ellipse, the moment as a function of 
angle is then extracted and decomposed into the Fourier 
series 



K(ip) = A Q + Ai sin(V>) + Bx cos(ip) 
+ A 2 sin(2V>) + B 2 cos(2V>) + . . . , 



(1) 



where the radial dependence of all A n 's and B n 's is im- 
plicit, since the above expression is for a single kinemetry 
ellipse. Here, ip is the azimuthal angle in the plane of 
the galaxy, measured from the major axis; points along 
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Fig. 2. — Left: Velocity field of a toy disk model (see Section 3.2), with kinemetry ellipses overlaid. One ellipse is emphasized with the 
solid line; the sampling of the rest of the ellipses are shown with the black dots. Center: The kinemetry expansion as a function of angle 
if) along the solid ellipse. The top panel shows the measured velocities (black points) and the fit with the B\ coefficient (red); the bottom 
panel shows the residuals from this fit (black points), and the higher order coefficients measured as a function of ip (A\ in red, A2 green, 
B2 blue, A3 magenta, and S3 cyan). Right: The kinemetry expansion from the center panel, now shown for all ellipses as a function of 
semi-major axis length a. The top panel shows Bi as a function of a; this is the rotation curve. The bottom panel shows the strength of 
the higher order coefficients (same colors as above), all of which are negligible, as would be expected in an ideal disk. In this toy model, 
deviations from zero in these coefficients reflect the noise in the velocity field. 



the ellipse are sampled uniformly in ip and are therefore 
equidistant if the ellipse is projected on to a circle. The 
series can be presented as a function of semi-major axis 
length a and in a more compact way, 



K(a,ip) = A (a) 



N 

E 



k n (a) cos[n(ip - <j> n {r)% (2) 



with the amplitude and phase coefficients (fc„,</> n ) de- 
fined as 



A n 



k„ = yjAl + Bl and 



6„ = arctan ( — ) . (3) 



The full moment map (i.e. velocity or velocity disper- 
sion) can thus be described by the geometry of the rings 
and the amplitude of the coefficients k n (or equivalcntly, 
A n and B n ) of the Fourier terms as a function of semi- 
major axis length a (Figure 2). 

To be more specific, the velocity field in an ideal rotat- 
ing disk is expected to be dominated by the cos ip term, 
since the velocity peaks at the galaxy major axis (ip = 0) 
and goes to zero along the minor axis (ip = Figure 
2). The power in the B\ term therefore represents the 
circular velocity at each ring a, while power in the other 
coefficients (normalized to the rotation curve, B\) repre- 
sents deviations from circular motion. In local galaxies 
with very high S/N observations, for example, various 
coefficients have been shown to identify bars/radial in- 
flow, lopsidedness/ warps, multiple components, and spi- 
ral structure (Schoenmakers et al. 1997; Wong et al. 2004; 
Krajnovic et al. 2006). In lower S/N data with sparser 
spatial sampling, these higher-order coefficients will also 
be affected by the rapid variations along each ring in- 
duced by the noise. 

The velocity dispersion field, on the other hand, is an 
even moment of the velocity distribution and, as such, 
its kinemetric analysis is identical to traditional surface 
photometry. In an ideal rotating disk, the velocity dis- 
persion will be constant along each ring and will decrease 
between rings of increasing semi-major axis length (Fig- 
ure 1). For this moment of the velocity distribution, 



the power in the Aq term as a function of semi-major 
axis length a will represent the velocity dispersion profile, 
and all azimuthally varying terms (higher-order kineme- 
try coefficients) will be zero. In analogy to the case with 
the velocity field, non-zero A n 's and i? n 's can thus iden- 
tify expected deviations from symmetry, in the form of 
lopsidedness and boxy/disky iso- velocity dispersion con- 
tours, but are also susceptible to variations caused by 
noise. 

2.2. Kinemetry of High-Redshift Systems 

Kinemetry was originally designed by Krajnovic et al. 
(2006) for use with very high S/N (> 100) stellar kine- 
matic data, as found in observations of bright galaxies 
in the local Universe. It has also been used to ana- 
lyze simulated mergers and merger remnants by Jesseit 
et al. (2007) and Kronberger et al. (2007). To apply this 
method to the much lower S/N emission line data ob- 
tained at z ~ 2, the breadth of the analysis must be 
somewhat restricted. We therefore employ kinemetry in 
a more limited capacity; rather than using kinemetry to 
measure and interpret subtle dynamical features of a ve- 
locity field, we instead use it to determine the strength of 
deviations of the observed velocity and dispersion fields 
from the ideal rotating disk case. This is identical to 
assuming that any deviations from the ideal case that 
might occur in a disk (e.g. lopsidedness, warps, spiral 
structure) induce less power in the higher Fourier coeffi- 
cients than those caused by the noise and much less than 
those that occur in a disturbed, merging system (com- 
pare Figures 2 and 3). 

The first step in the analysis is locating the center of 
the system, around which the kinemetry ellipses are con- 
structed. The robustness of this step is critical to the re- 
sult, since Krajnovic et al. (2006) show that an incorrect 
assumed center induces artificial power in the derived 
kinemetry coefficients. The primary result of a miscen- 
tering are elevated A2 and B2, but other coefficients (Aq, 
Ai, A3, and B3) are affected as well. It is therefore im- 
portant that we use a robust definition of the system 
center, such that the center of an ideal disk is accurately 
recovered. 
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This is a non-trivial task in the clumpy and irregu- 
lar Ha intensity distribution, the spatial distribution of 
which corresponds only to regions with enhanced star for- 
mation rates, is influenced by extinction, and therefore 
does not necessarily reflect the intrinsic mass distribution 
in the system (Figure 1). It is also not straightforward to 
robustly derive the location of the galaxy center from the 
velocity and velocity dispersion maps themselves (Kra- 
jnovic, private communication). We therefore take ad- 
vantage of the continuum distribution, which can be de- 
tected in the integral field data at a somewhat lower S /N 
underneath the emission lines. In general, the continuum 
surface brightness distribution is distinct from that of the 
emission line intensity (Figure 1 , see also examples in lo- 
cal spiral galaxies in Daigle et al. 2006) and has sufficient 
S/N to differentiate regions of strong continuum emission 
from those of weaker continuum emission. 

At z ~ 2, the detected continuum in near-infrared 
observations corresponds roughly to rest-frame i?-band. 
Observations at this wavelength therefore typically pro- 
vide an accurate probe of the stellar distribution, al- 
though the effects of extinction have been shown to be 
significant in some local dusty mergers (Arribas & Col- 
ina 2003). While it is difficult to know precisely how 
large a role extinction plays in the observed i?-band con- 
tinuum distributions of high-z observations, some con- 
straint can be provided from visual analysis of the data 
themselves. Several systems observed at z ~ 2 show a 
compact bright, central region of continuum emission, 
which remains single-peaked even at < 0'.'5 (4 kpc) res- 
olution (e.g. Figure 1), suggesting that this is an unob- 
scured measure of the stellar (and mass) distribution in 
these systems. 

As such, the continuum surface brightness distribution 
can be used to locate the center of a system. The natural 
choice for a galaxy's center is the position of the peak 
in the continuum intensity, which in an ideal disk also 
corresponds to the centers of the kinematic fields (Figure 
1); however, in noisy data, this must be carefully defined. 
We identify the center of the continuum distribution by 
using the data in the brightest 25% of the pixels and 
finding the continuum-intensity- weighted average of their 
positions. This definition of the system center thus has 
the added benefit that the center of an early-stage major 
merging system (with two distinct nuclei) will be directly 
between the two components and not skewed towards 
only one of the two mass concentrations. 

In high S/N data, the location of the center is the only 
parameter that must be determined prior to running a 
kinemetric analysis; the relevant geometric parameters 
(position angle and inclination) are derived during the 
kinemetry analysis. However, the lower S/N and coarse 
spatial sampling of the z ~ 2 data (see Figure 1) renders 
kinemetry's radius-by-radius solution for position angle 
and inclination rather unstable. For these systems, a 
much more robust solution for these parameters is found 
by considering the entire velocity field at once and solving 
for a global position angle and inclination. In nearby spi- 
ral galaxies, these quantities are observed to vary slowly 
throughout the system (Wong et al. 2004), thus making 
global values for position angle and inclination decent 
approximations. 

To determine the position angle and inclination of a 
system, we use the known effects of errors in these pa- 



rameters on the kincmctry coefficients A n and B n . Kra- 
jnovic et al. (2006) demonstrate that, in the kinemetric 
expansion of a velocity field, a slightly incorrect assumed 
position angle generates excess power in the A\, A 3} and 
i? 3 coefficients, while a slightly incorrect assumed incli- 
nation affects primarily the B$ term. However, large 
errors in these values can produce significant power in 
other coefficients as well. We therefore use the measured 
power in all coefficients to derive these parameters. 

We solve first for the global position angle by step- 
ping through all possible values, in increments of 3°, 
and performing a kinemetric expansion of the velocity 
field at each assumed value. Since this procedure fo- 
cuses on locating the angle of steepest velocity gradient, 
the axial ratio of the ellipses (i.e. the inclination of the 
system) does not affect the results and is therefore held 
constant at unity, forcing the kincmctry ellipses to be cir- 
cles. The goodness-of-fit of the assumed position angle is 
determined as the sum of the squared residuals between 
the "circular" velocity field (the 2D image reconstructed 
from only the B\ term) and the observed velocity field. 
These residuals reflect the combined powers in the higher 
coefficients, which appear as asymmetries in the velocity 
field. The curve of the goodness-of-fit as a function of 
position angle is then smoothed by 3 data points (=9°) 
to eliminate spurious results induced by the noise in the 
data, and the best-fit position angle is identified as the 
position angle that minimizes the smoothed curve. 

Assuming this position angle, we then solve for the 
global inclination of the system, again using kinemet- 
ric expansion of the velocity field. We test ~ 50 values 
for the inclination, evenly spaced between axial ratios of 
0.1 and 1.0. At every assumed ellipticity, kinemetry is 
performed, with the position angle held constant at the 
previously-determined best-fit value and the inclination 
held constant at the assumed value. The goodness-of-fit 
is determined as above, using the residuals between the 
"circular" velocity field and the actual measured values. 
The minimization of this method therefore finds the in- 
clination that produces the smallest deviations from the 
circular term. However, since ellipses of similar axial 
ratios to that measured from the morphology will max- 
imize the kinemetry coverage of the velocity field, and 
thus minimize the residuals measured near the edges of 
the system, this process induces a slight bias towards an 
inclination similar to that of the morphology. In prac- 
tice, this bias partially mediates the adverse affects of 
beam smearing, in which the circular beam increases the 
opening angle of the isovelocity contours of the inclined 
velocity field. As with the position angle determination, 
we smooth the goodness-of-fit curve by 3 data points 
to reduced the impact of the noise in the velocity field, 
and we find the inclination that minimizes the smoothed 
residuals curve. 

With the position angle and inclination held constant 
at these best-fit values, we can then perform a robust 
kinemetric analysis of the low-S/N high-z velocity and 
velocity dispersion fields. As in Krajnovic et al. (2006), 
kinemetry is carried out to the fifth order terms, and 
£>5, which include most of the physical deviations from 
symmetry in the field but are not overly affected by rapid 
point-to-point variations induced by noise. Since the 
kinemetry coefficients are by definition orthogonal, this 
somewhat arbitrary choice of where to stop the expan- 
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Fig. 3. — Same as Figure 2, but for an observed merger (IRAS 15206+3342) simulated at z ~ 2 (see Section 3.1). The irregular velocity 
field (left) is not well fit by a cos-0 term and therefore has significant power in all higher order coefficients (center). The combined effects of 
the lack of power in B\ from the poor fit to cos-0 and the power in all other coefficients (right) will produce normalized coefficients k n /Bi 
that are much higher than those of a disk (Figure 2). 



sion will not alter the measured values of the kinemetry 
coefficients. 

2.3. Criteria for Differentiating Disks and Mergers 

The main component of the analysis thus consists 
of running kinemetry, with pre-defined system center, 
global position angle, and global inclination, on both the 
velocity and velocity dispersion field of a galaxy. Since 
the quantity of interest is the deviations of these fields 
from those of ideal disks, we find the average deviation, 
defined as the average of the kinemetry coefficients that 
would be identically zero in a noiseless ideal disk. 

For a velocity field in an ideal, rotating disk, the only 
non-zero kinemetry coefficient ought to be B\, which we 
denote here as B\ v to indicate that this is the kineme- 
try result from the velocity field (Figure 2). Information 
about additional asymmetries are therefore contained in 
the higher-order terms, k 2tV ~~ fcs.u- Although Krajnovic 
et al. (2006) use only odd kinemetry terms (e.g. k 5 ) 
to describe odd moments of the velocity distribution 
(velocity field), we include both the even and the odd 
terms in the kinemetry expansion here, since mergers 
produce extremely disturbed velocity fields, with power 
in all kinemetry coefficients (Figure 3). In principle, one 
could also include the A\ )V term as a measure of asym- 
metry, since this term represents any velocity gradients 
(inflows/outflows) found along the minor axis. In prac- 
tice, however, we choose not to use the presence of a non- 
zero A% tV coefficient in the definition of non-ideal disks, 
since significant radial flows may be the result of such 
phenomena as outflows related to AGN/starburst winds 
or inflows induced by bar instabilities (see discussion in 
Section 5.2), and so do not provide information as to 
whether a system is undergoing a merger. We therefore 
compute the average k avgtV = (k 2 ^ v + + k 4:V + k 5yV )/A 
as a measure of the non-ideal components of a system. 

This average deviation k avg ^ v is normalized to the ro- 
tation curve B\ :V in order to assess the relative level of 
deviation, following the prescription of Krajnovic et al. 
(2006). This normalization has the added benefit of ac- 
counting for the loss of dynamic range in k avg ^ v at high 
redshifts, where the dim outer regions of a galaxy are too 
faint to be detected. In merging systems, greater radial 
coverage will correspond to greater detectable deviations 
from the ideal disk geometry, and it will also correspond 



to a greater dynamic range in the velocity gradient B\ 
(Figure 4). By normalizing k avg ^ v to Bx >v , then, a sys- 
tem will roughly retain the same amount of detectable 
asymmetry regardless of the radial extent of the data. 
We therefore define the asymmetry, or level of deviation 
from an ideal disk, in the velocity field to be 

_ / k a yg V \ . 

Vasym \ — ^ / > (^) 
\ B l,v I r 

where the average is over all radii (relative to the contin- 
uum peak) of the kinemetry ellipses. 

For the velocity dispersion field, the only non-zero 
kinemetry coefficient in an ideal, rotating disk is Aq^, 
which quantifies the velocity dispersion profile. All 
higher-order coefficients k\^ a - k^^ a therefore measure 
any asymmetries in the field. For this moment of the 
velocity distribution, then, k avg ^ = (k\^ a + k 2 , a + ^3,0- + 
fc+cr + fc5,o-)/5 will contain information about the devia- 
tions from the ideal case. 

Krajnovic et al. (2006) note that, when looking at the 
velocity dispersion of the stellar component, the appro- 
priate normalization for this even moment is the Aq <(T 
coefficient. In the stellar case, the Ao >tT term traces the 
mass of the system; however, this is not generally true 
of the velocity dispersion of a gas component, which can 
also be affected by shocks, especially in the violent envi- 
ronment of a major merger. In the kinemetry of gas 
kinematics, if there is significant rotation, the poten- 
tial is often more reliably probed by the rotation curve, 
B\. v . Furthermore, in Figure 4, the weakness of using a 
normalization to Aq^ in the presence of varying radial 
coverage is clear. As with the velocity field deviations 
k avgiV , the velocity dispersion deviations k avg ^ become 
stronger with radius; a loss in radial coverage therefore 
directly corresponds to a loss in dynamic range in the 
asymmetries. In contrast, the value of the Aq i(t circu- 
lar term remains roughly constant, even when observed 
with the broader PSF of z ~ 2 observations. The value 
of kayg^o- 1 A 0t(7 consequently decreases significantly for a 
given system when observations are less sensitive to the 
outer regions. A more appropriate normalization is the 
rotation curve of the velocity field, B\^ v , which is both a 
more reliable measure of the system's mass and respon- 
sive to the loss of dynamic range with decreased sensi- 
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Fig. 4. — From left to right: The map of a given velocity moment, followed by the circular term of its kinemetry expansion, the higher 
order coefficients of the expansion k avg , and the derived asymmetry parameter, all as functions of radius from the system center. This 
is shown, from top to bottom, for the velocity and velocity dispersion of the z ~ 0.1 ULIRG Mrk 273, and for the velocity and velocity 
dispersion of the simulated version of this system at z ~ 2 (see text for details). In the velocity maps, the diminished spatial extent of 
the high-redshift data corresponds to less resolution of both the rotation curve and the higher order coefficients. On the other hand, the 
velocity dispersion maps show negligible change in the dispersion profile with changing spatial coverage, while the higher order terms, as 
with those from the velocity field, lose dynamic range. 



tivity. We therefore define the asymmetry in the velocity 
dispersion field as 



1 asym 



(5) 



where the average over all radii is unaffected by the com- 
bination of velocity and velocity dispersion kinemetric 
coefficients, since the kinemetry ellipses for the two maps 
are identical by the construction described in Section 2.2. 

3. APPLICATION TO TEMPLATE GALAXIES 

To assess the capabilities of these criteria, we draw on a 
sample of observed low-redshift disks and mergers, which 
are then "observed" as if they were at z ~ 2, as well as 
on a sample of synthetic high-z systems, created from 
simulations of varying complexity. This set of template 
galaxies spans a large range in morphology, merging his- 
tory, nuclear activity, and star-formation rate. 

3.1. Observed Systems 

Our template local galaxies are drawn from two sam- 
ples: the SINGS spiral galaxy survey (z ~ 0; Kenni- 
cutt et al. 2003), as observed in Ha by Hernandez et al. 
(2005); Daigle et al. (2006); Chemin et al. (2006), and the 
Ha observations of low- z ultraluminous infrared galaxies 
(ULIRGs, z ~ 0.1) of Colina et al. (2005). Together, 



these data sets provide high-quality local observations of 
disks and mergers, respectively. 

The SINGS galaxies were observed in Ha emission 
with the FANTOMM Fabry-Perot scanning interferome- 
ter (Hernandez et al. 2003). The template galaxies used 
here were all observed at the 1.6m telescope of the Obser- 
vatoire du mont Megantic in Quebec, Canada (Hernan- 
dez et al. 2005; Daigle et al. 2006; Chemin et al. 2006). 
The pixel scale of the velocity maps is 1'.'61, and the spec- 
tral resolution i?=12200-20000. From this sample, we 
use six spiral galaxies (NGC 925, NGC 3198, NGC 4321, 
NGC 4579, NGC 4725, and NGC 7331), which together 
span a range of Hubble types, nuclear activity, star for- 
mation rates, and distribution of star-forming regions 
throughout the system. 

The ULIRG merger systems were observed in Ha emis- 
sion and optical continuum with the INTEGRAL fiber- 
fed integral field system (Arribas et al. 1998), mounted on 
the 4.2m William Hershel Telescope of the Roque de Los 
Muchachos Observatory of La Palma, Spain. These data 
have resolution elements (fibers) of 0'.'9 diameter and 
spectral resolution i?=1500 (Colina et al. 2005). From 
this sample, we use eight mergers (Mrk 273, Arp 220, 
IRAS 08572+3915, IRAS 12112+0305, IRAS 14348- 
1447, IRAS 15206+3342, IRAS 15250+3609, and 
IRAS 17208-0014), which, in analogy with the spiral sam- 
ple, span a range of merger stage, nuclear activity, star 
formation rates, and distribution of star-forming regions 
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throughout the system. 

For both the disk and merger template galaxies, we 
simulate observations of these systems as if they were 
at z ~ 2 and were observed with the SINFONI integral 
field unit (Eisenhauer et al. 2003; Bonnet et al. 2004), 
on the Very Large Telescope (VLT) at Cerro Paranal, 
Chile. To artificially "redshift" the template galaxies, 
we convolve the datacubes to the mean seeing of such 
observations (FWHM of 0'.'5 ~ 4 kpc; Forster Schreiber 
et al. 2006 achieved this beam size as typical seeing, while 
Genzel et al. 2006 used adaptive optics to reach a see- 
ing FWHM of 0'.'15) and to the spectral resolution of the 
instrument (75 km s _1 at 2.2/im), interpolate the data 
onto the larger pixel scale (spatial scale of 0'.'125 ~ 1 kpc; 
velocity scale of 33 km s _1 ), account for the cosmologi- 
cal surface brightness dimming, and add Gaussian noise 
such that the S/N of the resulting datacubes is compa- 
rable to that of VLT/SINFONI observations. The Ha 
kinematics are extracted from these datacubes using the 
same technique as that of Forster Schreiber et al. (2006) 
and Bouche et al. (2007) on actual SINFONI observa- 
tions and are shown in Appendix Figures 8 and 9. 

The resulting sample of observed galaxies thus totals 
14; six of local disks and and eight of local mergers, all 
"redshifted" to z ~ 2. 

3.2. Simulated Systems 

Because this naive "rcdshifting" of local systems can- 
not truly emulate the proto-galaxy population at z ~ 2, 
in which such phenomena as massive (10 8 ~ 9 M©) clumps 
of star formation affect the structure and kinematics of 
the systems (Bournaud et al. 2007), we supplement our 
template galaxies with synthetic systems of various com- 
plexity. For this, we use both toy disk models, in which 
we have a complete understanding of all aspects of the 
data, as well as the detailed hydrodynamic cosmologi- 
cal simulations of Naab et al. (2007), from which we can 
"observe" synthetic systems at z ~ 2. 

For the toy disk models, we use the modeling routines 
described in Forster Schreiber et al. (2006), which gen- 
erate simple models of azimuthally symmetric rotating 
disks, parametrized by mass, inclination, scale length, 
scale height, and isotropic velocity dispersion. Appro- 
priate pixel sizes, beam smearing, spectral resolution, 
and noise levels are also included in the construction 
of the toy datacube. We improve on this model by al- 
lowing both an azimuthally-symmctric coupled mass and 
light distribution, and an overlaid light-only distribution 
(not azimuthally symmetric). The former therefore rep- 
resents the stellar component of the system, while the 
latter produces very bright regions corresponding to ar- 
eas of increased star formation. This addition allows us 
to generate models that more realistically simulate the 
clumpy and irregular Ha intensity distribution of systems 
at z ~ 2. We create five models with this technique, 
all of which have a centrally-peaked and azimuthally- 
symmetric mass distribution. We vary their light distri- 
bution (star-forming regions) as follows: one galaxy is 
nearly edge-on and is azimuthally symmetric around a 
centrally-peaked light distribution, one is nearly edge-on 
with a varying light distribution, one is of intermedi- 
ate inclination and is azimuthally symmetric around a 
centrally-peaked light distribution (shown in Figure 1), 
one is of intermediate inclination with a light distribution 



much more extended on one side of the galaxy, and the 
last is of intermediate inclination with the light distribu- 
tion illuminating only a single side of the galaxy. From 
these models, we extract the emission line kinematics and 
continuum intensity by fitting a Gaussian with a constant 
offset to each spectrum (see Appendix Figure 10 for the 
models and their derived kinematics). 

Our cosmologically-simulated models are drawn from 
Naab et al. (2007). The synthetic z ~ 2 systems come 
from high resolution (10 6 gas particles and 10 6 dark mat- 
ter particles per halo) SPH simulations of 8 halos with a 
variety of mass accretion histories. Based on these mass 
accretion histories, we select time snapshots in which the 
halos have been evolving without a major merger in re- 
cent history or in the near future as well as snapshots 
in which the halos are visibly in the process of a major 
merger. For the latter group, this selection criteria es- 
sentially require the system to have a double nucleus, 
separated by < 5 kpc. From the 8 halos, we find 6 
snapshots between z = 1.8 and z — 2.8 in which the 
systems are unambiguous disks and 4 snapshots between 
z = 1.6 and z = 3.0 in which the systems are unambigu- 
ous mergers. These 10 snapshots of the model galax- 
ies are "observed" by converting the star-formation rate 
to Ha emission using the conversion factor from Kcnni- 
cutt (1998), accounting for cosmological surface bright- 
ness dimming, binning to SINFONI-size pixels, convolv- 
ing the data to the appropriate spatial and spectral res- 
olutions, and adjusting the (Gaussian) noise level such 
that the S/N in Ha is comparable to that of the SINS 
observations. The continuum intensity for each system 
is "observed" through a similar process, using the stel- 
lar mass and converting to i?-band luminosity with an 
assumed M/Lr — 1, typical for a Kroupa IMF in star- 
forming systems at this redshift (Fontana et al. 2004). 
We then rotate the halo to a random inclination, as given 
by a sin(i) probability distribution function, and extract 
the Ha kinematics from the resulting datacube (see Ap- 
pendix Figures 11 and 12). 

The resulting sample of model galaxies totals 15; five 
toy disks, six cosmologically simulated disks at z = 1.8 — 
2.8, and four cosmologically simulated mergers at z = 
1.6-3.0. 

3.3. Classification 

We test the criteria described in Section 2.3 on the 29 
template galaxies to determine how well we can differen- 
tiate disks and mergers based on warm gas kinematics. 
For each system, we perform the analysis of Sections 2.2 
and 2.3 to measure v asym and a aS ym for these templates. 
(The results of the kinemctric analysis are summarized 
in Appendix Table 2 and shown in Appendix Figures 8 
through 12.) Since this method does not lend itself to 
straightforward error propagation, we use Monte Carlo 
realizations of the noise in the data to measure the prob- 
ability distribution functions (PDFs) of the asymmetries 
in these systems. 

For each template system, the Monte Carlo realizations 
consist of creating 1000 different realizations of the mo- 
ment maps - the continuum intensity, and the emission 
line intensity, velocity, and velocity dispersion - based on 
their corresponding error maps. These error maps cor- 
respond to the measurement errors of the velocity mo- 
ments, as derived when fitting the kinematics from the 
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Fig. 5. — Asymmetry measure of the velocity and velocity dispersion fields for (left) all of the template galaxies and (right) all of the 
template galaxies except the obviously misclassified ULIRG, IRAS 12112+0305. The probability distributions in this space are shown with 
shading for the template disks (blue) and mergers (red), as derived from the Monte Carlo realizations. Inset are the PDFs for the total 
kinematic asymmetry (K asyrn ) for disk and mergers. The empirical delineation of K asyrn = 0.5 cleanly separates the two classes, as is 
especially visible in the right panel. 



asyrn 



a asym = 0.5, as visible 



datacubes. For each moment map, we perturb the ob- 
served data points by randomizing them, using Gaussian 
noise parametrized by the measured (1-cr) errors. The 
new set of maps is then rerun through the entire analysis 
described in Sections 2.2 and 2.3. 

Fi gure 5 illustrates the resulting v aS y m and <J a sym mea- 
surements for the template systems. In this figure, all of 
the results from the Monte Carlo realizations are plot- 
ted, with red shading indicating the resulting PDF of 
the merger templates and blue shading indicating that 
of the disk templates. These two classes can be cleanly 
separated by the empirical delineation of total kinematic 

asymmetry K asym = 

in the inset. 

The majority of the disks (89%) are located in the 
lower left (low v aS ym, low a asym ) of the diagram, with the 
small deviations from the ideal case (v asyrn = 0, a asym 
= 0) coming from noise, thickness of the disk, and other 
kinematic features such as warps and multiple compo- 
nents. The mergers, for the most part, show strong de- 
viations from zero in both v asym and a asym and are lo- 
cated towards the upper right of the plot. However, 20% 
of merging systems remain indistinguishable from disks. 
This is largely due to a single ULIRG, IRAS 12112+0305, 
whose velocity and velocity dispersion fields appear reg- 
ular at z ~ 2 (see Appendix Figure 9), although other 
systems contribute to a lesser extent as well (Figure 5). 
Based on these results, we can roughly estimate the er- 
rors in these criteria and can expect to correctly classify 
~ 89% of disks and ~ 80% of mergers. 

Because these conclusions are based on a detailed and 
complicated analysis (Section 2.2), we tested on several 
systems how changes in the assumed center, position an- 
gle, and inclination would affect the classification of the 
system as a disk or a merger. We first examined varia- 
tions in the assigned center and found that, for disks, the 



classification of a system as such is virtually independent 
of the center, except in extreme miscenterings when the 
center is assigned to the very edge of the system. How- 
ever, the location of the center is more important in the 
case of mergers with a double nucleus. In these systems, 
if the center is skewed too far towards one of the two 
mass concentrations (both of which have some ordered 
rotational motion), the system can be misclassified as a 
disk. This reinforces the necessity of choosing the cen- 
ter of mass of the system, via the continuum-intensity- 
weighted center as described in Section 2.2. 

The test of variations in position angle indicated that 
the classification of a system as a disk or merger is even 
more robust against changes in this parameter. In the 
case of disks, these systems maintain their low K asym 
through extreme variations in position angle, until the 
position angle is aligned within ~ 10° of the minor axis. 
In mergers, the position angle has no physical meaning 
but, by the process described in Section 2.2, has been 
defined such that it minimizes K asym . As a results, in 
this systems, all variations in position angle only increase 
their measured asymmetries, making their classification 
as mergers even stronger. Finally, variations in inclina- 
tion were also examined and were shown in nearly all 
cases to have no effects on the classification of either 
disks or mergers. These tests of the effects of assumed 
inclination, position angle, and center on the classifica- 
tion of a system therefore suggest that, even in systems 
where there remains some uncertainty about the values 
of these parameters, our method will provide a robust 
classification with K asym . 

4. APPLICATION TO HIGH-REDSHIFT GALAXIES 

Having tested these criteria on our template galaxies, 
we now apply them to the 'unknowns,' the high-z sys- 
tems observed in the SINS program (Forster Schreiber 
et al. 2006; Genzel et al. 2006; Bouche et al. 2007; Cresci 
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Fig. 6. — Left: Comparison of the "redshifted" template galaxies' typical errors in velocity and velocity dispersion (dashed lines) with 
the measured median errors of each of the SINS program galaxies. SINS systems with smaller errors than those of the templates are shown 
in red. Right: The horizontal axis compares the median number of beams covering each system in the template galaxies (vertical line) with 
those of the SINS galaxies. The vertical axis compares the median number of beams covering data with continuum S/N> 3, as used to 
determine the galaxy centers, in the template galaxies (horizontal line) with those of the SINS systems. Red points are the same galaxies 
shown in red in the left panel; red and blue points are those systems with data of better quality than that of the templates, based on both 
panels. Only these systems are included in our analysis, with one additional inclusion. We also analyze BzK-15504, shown with the red 
and green point, which fulfills three of the four data quality requirements but has insufficient S/N in its continuum emission. The lower 
S/N per pixel is in part due to the finer spatial sampling of this system, which Genzel et al. (2006) observed with the O'.'l pixel scale. 
Fortunately, the continuum emission of this system has been observed in deep VLT/NACO imaging with laser guide star adaptive optics, 
which revealed a distribution consistent with that measured from the SINFONI data. We are therefore confident that the continuum center 
measured from the SINFONI data does in fact reflect the center of the mass distribution in this system. 



et al. in prep). Here we examine those galaxies with suf- 
ficient spatial resolution and data quality to perform the 
kinemetric analysis described above. Galaxies that are 
unresolved, have large errors in the kinematics, or have 
significantly lower S/N than the majority of the sample 
are omitted from the present discussion; in future work, 
we will expand our criteria to include variations in these 
quantities. 

4.1. Data 

The SINS z ~ 2 sample is taken from large 
photometric samples, in which high-z objects are 
identified through either their rest-frame ultraviolet 
color/magnitude (BM/BX criterion: Adelberger et al. 
2004; Steidel et al. 2004; Erb et al. 2006a,b), or their rest- 
frame optical properties (s-BzK: Daddi et al. 2004a, b; 
Kong et al. 2006, GDDS: Abraham et al. 2004). These se- 
lection criteria sample luminous (L ~ 10 n ~ 12 L©) galax- 
ies with an range of star formation rates (SFR ~ 10 — 200 
M /yr) and ages (50 Myr - 2Gyr; Erb et al. 2006a; Daddi 
et al. 2004a, b). From these photometric samples, our se- 
lection criteria emphasize somewhat brighter (<F(Ha)> 
of 10~ 16 compared to 6 x 10 - erg/s/cm 2 ) systems, 
with broader line widths (<v c > of 175 ± 68 compared 
to 140 km/s) than the average galaxy in the Erb et al. 
(2006b) sample. Both samples have similar mean dy- 
namical masses. 

The SINS galaxies were observed in Ha emission, 
which at z ~ 2 is redshifted to the if-band, with 
VLT/SINFONI. Approximately i?-band continuum emis- 
sion is visible beneath the strong emission lines. Most 
of the data have 0'.'125 x 0'.'250 pixels, sampling a typ- 
ical PSF FWHM of 0'.'5, and have a spectral resolution 



of R ~ 4000. Additionally, a few systems have been 
observed with adaptive optics, enabling the use of the 
finer pixel scale of 0'.'05 x O'.'IO to sample the typical PSF 
FWHM of 0'.'15. The data analysis is described in Forster 
Schreiber et al. (2006); here we add the additional step of 
spatially binning the datacube to a minimum amplitudc- 
to-noise ratio of 5 with the Voronoi binning technique of 
Cappellari & Copin (2003). This reduces the spatial res- 
olution at the fainter edges of a system, where several 
spatial elements must be summed, but the amplitude-to- 
noise requirement ensures a more robust measurement of 
the kinematics in each bin. 

The SINS systems were observed with integration 
times ranging from 1.5 hours to over 8 hours (Forster 
Schreiber et al. 2006; Genzel et al. 2006; Bouche et al. 
2007; Cresci et al. in prep), which, when coupled with 
differences between systems, results in a wide range of 
errors on the kinematic measurements, S/N levels, and 
number of beams covering the systems. In Figure 6, we 
compare the assumed values for our template galaxies 
(typical values for the entire sample) to the actual val- 
ues for individual SINS galaxies (see figure caption for 
details). For this sample, the most critical and most lim- 
iting requirement is that the systems are well resolved, 
which we define as covered by > 3 beams, although we 
also require that the continuum be detected with a sig- 
nificant S/N level and that the kinematic measurements 
be sufficiently precise. Analyzing data that do not meet 
these standards will adversely affect the spatial sam- 
pling of the kinemetry ellipses, their centering, and the 
precision of the derived asymmetry measures v asym and 
<j aS ym.'i the calibration of our analysis to such lower qual- 
ity data is outside the scope of this paper. We find that 
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5.3 


optical 


D3a-6397 


1.51 


65 


6.3 


7.6 


optical 


K20-ID7 


2.22 


84 


3.5 


4.5 


optical 


K20-ID8 


2.22 


42 


1.8 


5.4 


optical 


Average 


2.12 


81 


6.8 


5.9 




SINS UV Sample 


2.26 


60±36 


4.1±3.7 


4.3±2.1 


UV 


SINS Optical Samplc f 


2.02 


46±49 


4.2±4.5 


4.7±1.6 


optical 



a As derived from the total Ha flux measured from the SINS data. These values are converted to SFR 
via the Kennicutt (1998) conversion, with an additional conversion to a Chabricr IMF. Extinction 
Ay is measured from SED modeling (sec Forster Schreiber et al. in prep), and the Calzetti extinction 
law is used to estimate extinction Ah a, from which the extinction correction is calculated. 

b From the measured circular velocity, R1/2, and axial ratio of the intensity distribution. 

c Rest-framc waveband in which the object was identified. 

d Observed with LGS, yielding a PSF FWHM of 0'.'15; all other systems were obtained in seeing- 
limited mode with a mean PSF FWHM of 0'.'5. 

e Systems observed in the SINS program that were originally identified with rest- frame UV pho- 
tometry. Values given are mean and standard deviation of the sample. 

f Systems observed in the SINS program that were originally identified with rest-frame optical 
photometry. Half of this sample has z ~ 2.3, and half 2 ~ 1.5. Values given are the mean and 
standard deviation of the sample. 



11 of the z ~ 2 systems have equal or higher quality data 
than the templates and thus their kinematic asymmetries 
can be reliably determined. We therefore perform our 
analysis on only these systems. 

The properties of this subsample are listed in Table 
1 and compared to the average properties of the UV- 
selected and optically-selected samples observed in the 
SINS program. In most cases, the star formation rates, 
dynamical masses, and half-light radii of the objects an- 
alyzed here are slightly higher than the mean values for 
these parameters from the full SINS sample but still 
within the standard deviations. Since the SINS sample 
contains a number of spatially unresolved compact ob- 
jects (Forster Schreiber et al. 2006), these statistics are 
consistent with the data requirements outlined above, in 
which only the well resolved (and thus large half-light ra- 
dius and likely high mass) objects are selected for analy- 
sis. It is these resolved objects in which the most progress 
can be made; for this reason, these systems were also the 
focus of Forster Schreiber et al. (2006) and Genzel et al. 
(2006). 

In their analyses, Forster Schreiber et al. (2006), Gen- 
zel et al. (2006), and Bouche et al. (2007) study the SINS 
data using the framework of rotating disks. These au- 
thors have argued from qualitative examination of the 
data, along with analysis of the rotation curves and 
comparison to model disks in several cases, that many 
of these UV/optical-selected systems are candidate disk 
galaxies, although they note that a minority do in fact 
have the disturbed velocity fields expected of mergers. 
We are now in a position to quantitatively test these 
conclusions. 

4.2. Classification 



We perform the analysis of Sections 2.2 and 2.3 on the 
SINS systems, with the same 1000 Monte Carlo realiza- 
tions, as done on the template z ~ 2 systems (Section 
3.3). In Figure 7 we plot the resulting v asym and a asym 
measurements for these galaxies. These results confirm 
the analyses of Forster Schreiber ct al. (2006) and Gen- 
zel et al. (2006) in that > 50% (8/11) of galaxies in the 
SINS subsample studied here are consistent with a ro- 
tating disk interpretation (see Appendix Table 2 for a 
list of results for each system). Given the error rate 
of these criteria, as found with the template systems, 
this suggests that these results sample a parent popula- 
tion composed of 8.4 disks and 2.6 mergers, such that 
89% x 8.4 + 20% x 2.6 = 8 systems are observed as disks 
and 80% x 2.6 + 11% x 8.4 = 3 systems are observed as 
mergers, making the fraction of disks in the subsample 
of Table 1 as high as 8.4/11 = 75%. Additionally, we 
expect that 20% x 2.6 ~ 1 merger will be misclassified as 
a disk and that 11% x 8.4 ~ 1 disk will be misclassified 
as a merger. 

It is of special interest to specifically investigate the 
most qualitatively convincing disks in the sample iden- 
tified by Forster Schreiber et al. (2006, SSA22a-MD41, 
Q2343-BX389, Q2343-BX610) and Genzel et al. (2006, 
BzK-15504). We highlight these systems in Figure 7 and 
indeed find that all four are consistent with a rotating 
disk interpretation. Forster Schreiber et al. (2006) also 
point out that two of their systems that may in fact be 
mergers; the one such system with sufficient data quality, 
Q1623-BX528, is identified in Figure 7 and is in fact a 
merger. Our criteria thus make intuitive sense based on 
visual analysis of the observed velocity and velocity dis- 
persion fields and therefore demonstrates quantitatively 
the validity of the rotating disk interpretation as applied 
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Fig. 7. — Asymmetry measure of the velocity and velocity dispersion fields for the SINS program galaxies that have high enough quality 
data for such analysis, overplotted on the disk and merger template PDFs from Figure 5. The line indicates the division between disks 
and mergers at K asyrn = 0.5. The probable disks identified by Forster Schreiber et al. (2006) and Genzel et al. (2006) are indicated here 
as green triangles. The merger identified by Forster Schreiber et al. (2006) is shown with the cyan square. Sample velocity fields of SINS 
disk-like and merger-like systems are shown at right; the full analysis of the SINS sample is recorded in Appendix Table 2 and shown in 
Appendix Figures 13 and 14. 



to the SINS sample. 

5. DISCUSSION 

Sections 3.3 and 4.2 illustrate the efficacy of our dy- 
namical criteria in identifying systems undergoing ma- 
jor mergers and their applicability to the SINS high-z 
survey; however, there are a number of caveats in this 
method that merit closer investigation. Here we expand 
the discussion of Section 4.2 to highlight the impact of 
this technique on high-z studies of galaxy formation, as 
well as address the most critical issues with this tech- 
nique: the ability of the Ha kinematics to reflect the 
underlying dynamics of a system, and the observational 
constraints on probing the structure of a distant system 
in detail. 

5.1. Implications for Galaxy Formation at 
High-Redshift 

Of the z ~ 2 systems analyzed here, we find that 8/11 
(> 50%) are consistent with a rotating disk interpre- 
tation. These results, coupled with the large stellar and 
dynamical masses of these systems {M dyn ~ 10 10 " 11 M ; 
Table 1), quantitatively confirm that massive disks were 
already in place at this redshift. It furthermore strength- 
ens the conclusions of Forster Schreiber et al. (2006) and 
Genzel et al. (2006) that the high (up to 100 M Q yr" 1 ) 
star formation rates of these systems are not the conse- 
quence of a recent violent merger but rather are happen- 
ing within disks dominated by ordered rotation. 



In the case of a typical massive (~ 10 11 M ) star- 
forming galaxy, BzK-15504, which Genzel et al. (2006) 
studied in detail, modeling of the spectral energy distri- 
bution (SED) from broad band photometry was used to 
measure the stellar mass and age of the system. With 
the stellar mass M* 8 x 10 10 M and the star formation 
rate of ~ 100 M yr -1 (as measured from Ha), together 
with the assumption of a constant star formation rate, all 
the stars in this system could have been formed in ~ 500 
Myr (Genzel et al. 2006). This number agrees well with 
the stellar age measured from the SED fitting (300 — 800 
Myr), making it likely that this system formed rapidly 
with continuous star formation (and therefore mass in- 
flow) at its current rate (Genzel et al. 2006). Given that 
the above analysis (Figure 7) indicates that this system 
has not undergone any recent major merger activity, and 
given that this system is typical for its population, our 
results provide new and direct empirical evidence that 
the smooth accretion mechanism can play an important 
role in the early stages of the evolution of massive galax- 
ies. 

The diagnostic tool described here is critical in expand- 
ing our understanding of structure formation and evolu- 
tion in the early Universe. With current extensive data 
sets of high-z systems, including broad-band photome- 
try and integral-field kinematic observations, much can 
be learned about the stellar populations, star formation 
processes, and nuclear activity of the galaxies evolving 
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in a critical epoch of the Universe's history. We now 
add another crucial piece to the study of high- z systems, 
by introducing a method to quantitatively evaluate the 
dynamical state of a system and therefore to link that 
system's observed properties with a major merger event 
or with a more quiescent evolutionary history. 

5.2. Ha as a Probe of a System's Dynamics 

The reliability of using Ha emission to study the struc- 
ture of a galaxy is, at first glance, rather unclear. The 
motions of the warm gas are not guaranteed to reflect 
those of the underlying stellar distribution, since the for- 
mer component is much more easily disturbed - with 
gravity or pressure fluctuations - than the more massive, 
collisionless stars. On the other hand, the dissipative gas 
component also more efficiently relaxes into a thin disk 
and could conceivably demonstrate ordered rotational 
motion while disturbances in the stellar (and mass) dis- 
tribution persist (e.g. Mihos & Hernquist 1996). These 
competing effects, together with the ambiguity caused 
by the limitations of spatial resolution, could conceiv- 
ably render it difficult to interpret the kinematics of the 
warm gas as uniquely representing a disk or a merger. 

Furthermore, the kinematics of the warm gas are also 
expected to reflect such phenomena as the large-scale 
gas flows that feed active nuclei and the powerful galac- 
tic winds from AGN/starburst activity, both of which at 
z ~ 2 arc thought to play an important role in regulating 
the star formation history of the Universe. For this rea- 
son, wc designed our criteria with input from templates 
that were likely to include as many of these phenomena 
as possible. Our sample of disks includes a kinemati- 
cally perturbed system, several barred systems, and two 
Seyfert galaxies. One of these active systems (NGC 4579) 
is a barred galaxy in which the radial motion of the gas 
is clearly visible as a strong velocity gradient along the 
minor axis (Gonzalez Delgado & Perez 1996; Daigle et al. 
2006). We also observe such an inflow in a SINS system, 
BzK-15504, as described by Genzel et al. (2006). In this 
system, a strong velocity gradient is also seen along the 
minor axis, in the form of a high A\. v term, presumably 
corresponding to the driving of fuel towards a growing 
bulge with an embedded active nucleus (Genzel et al. 
2006). We can nevertheless robustly identify both the 
"redshifted" NGC 4579 and BzK-15504 as disks, since 
the kinemetric signature of the inflows is restricted to 
the A\_ v coefficient, which is excluded from our analysis 
for precisely this reason (and could potentially be used in 
future work as a tracer of AGN feedback at high redshift). 
As shown in Section 3.3, this omission does not diminish 
our ability to detect mergers, since the disturbed veloc- 
ity fields of such systems also produce power in all the 
higher-order kinemetry coefficients. 

To discount the possibility of contamination by large- 
scale outflows in the z ~ 2 kinematics, we compare with 
observations of local systems with substantial galactic 
winds. In general, line emission from these "superwinds" 
in local starbursting galaxies dominates neither the total 
Ha surface brightness nor the emission line kinematics 
along the galaxy major axis (Lehnert & Hcckman 1995, 
1996a, b). The gas flowing outward along the minor axis 
has been shown to contribute only a small fraction of the 
total Ha luminosity; in even the extreme cases of local 
infrared-luminous and highly-extincted starburst galax- 



ies, Armus et al. (1990) show that the extended emission- 
line gas outside of a few kpc from the nucleus accounts 
for < 25% of the total Ha emission. In a few cases in 
our SINS sample, our data reveal a high- velocity outflow 
component, but this wind contributes < 10% to the to- 
tal Ha luminosity of the system. Similar results have 
been found with integral field data of a z ~ 2 SMG by 
Nesvadba et al. (2007). 

For the majority of the SINS z ~ 2 systems, then, it is 
likely that we can rule out possible contamination from 
superwinds and that we have sufficiently accounted for 
the affects of large-scale gas inflows. The z ~ 2 data 
are thus more consistent with the scenario in which the 
emission line kinematics trace dynamical features similar 
to those of the template galaxies. 

This argument suggests that the kinemetric analysis 
of a system is largely independent of SFR-driven phe- 
nomena and therefore of the overall SFR, a fact that is 
born out by the data themselves. In local spiral galaxies, 
where the relatively low SFRs (~ 1 — 10 M yr -1 ) re- 
flect low gas fractions, the dynamics of the warm gas are 
driven by those of the dominant, rotating stellar compo- 
nent. In contrast, the SINS galaxies, with SFRs charac- 
teristic of ULIRGs (up to and exceeding 100 M Q yr _1 ), 
have significantly higher gas fractions (f gas ~ 0.4; Tac- 
coni et al. 2006; Bouche et al. 2007) and consequently a 
dynamically important dissipative component. The gas 
in these systems quickly concentrates into massive, pow- 
erful star- forming regions, which drive turbulent motion 
and thus high velocity dispersions in the warm gas. How- 
ever, the clumps detected in SINS systems have not fa- 
tally disrupted the dynamics of their host galaxy; > 50% 
of SINS systems display velocity fields consistent with 
ordered rotation and regular (though elevated) velocity 
dispersions. Despite orders of magnitude difference in 
SFR between local spiral galaxies and SINS systems, it 
appears that the processes governing the Ha kinemat- 
ics in the two regimes are remarkably similar. This is 
consistent with results from Bouche et al. (2007), who 
demonstrated that star formation at high SFR and high 
redshift is governed by the same physics as in the lo- 
cal Universe. For the purpose of the analysis presented 
here, then, Ha emission can be effectively employed as a 
tracer of a system's dynamics for a large range of SFR 
and redshift. 

5.3. Continuum Surface Brightness Distribution 

One of the unique capabilities of the integral-field data 
used in our analysis, beside providing spatially resolved 
kinematics of a system, is separating the emission com- 
ing from star-forming regions, as visible in Ha emission, 
from that of the underlying stellar background, visible in 
approximately i?-band continuum emission. Morpholog- 
ical criteria using the continuum distribution to detect 
mergers have been implemented at lower redshift by e.g. 
Conselice (2003) and Lotz et al. (2004) using HST data. 
However, at z ~ 2 the stellar distribution is difficult to 
probe at optical (rest-frame UV and B) wavebands, since 
this emission can be severely affected by extinction in 
gas- rich systems (e.g. Colina et al. 2005). It is therefore 
of interest to examine the rest-frame i?-band continuum 
emission measured directly in our integral field data. 

Unfortunately, these methods require both high S/N 
and resolution elements smaller than 1 kpc/pixel, which 
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with SINFONI at z ~ 2 is attainable only in AO-assisted 
observations. Since our observations were optimized for 
the analysis of the line emission properties of the sources 
and were mostly obtained in seeing-limited mode, we 
tested several simpler criteria, based on kinemetry of 
the continuum distribution (which by definition is iden- 
tical to surface photometry). Any asymmetries in the 
continuum distribution would likely be due to the pres- 
ence of multiple mass concentrations, as in an early-stage 
merger. We find that the systems identified as potential 
early-stage mergers are, as expected, a subset of those 
identified as mergers by the kinematic criteria. This re- 
sult highlights the unique capabilities of directly probing 
the dynamical properties of high-z systems with integral 
field observations. 

We do not develop further any analysis of the con- 
tinuum emission in our SINS IFS data, as the low 
S/N of this emission limits the strength of the conclu- 
sions that can be drawn from such analyses. Forth- 
coming sensitive high-resolution (HST and ground-based 
AO-assisted) imaging of the SINS systems in the near- 
infrared will provide the necessary resolution and S/N 
required for the quantitative morphological analyses of 
broad-band emission (Forster Schreiber et al. in prep). 
The comparison of such results with those from the kine- 
metric analysis developed here will provide a valuable 
further probe into the nature of high- z systems. 

5.4. Limitations of the Method 

For the analysis described here, we designed our tem- 
plate systems to have S/N, resolutions, and spatial ex- 
tent comparable to typical values of the SINS survey at 
z ~ 2. In Section 4.1, we illustrated that a subsamplc 
of 11 SINS observations have higher quality data than 
the templates and thus have morphologies that can be 
reliably measured. 

For the remainder of the systems, and for systems 
at other redshifts or measured under different observing 
conditions, the current criteria cannot be blindly applied. 
We expect that there will be fundamental data limita- 
tions beyond which these criteria cannot distinguish disks 
from mergers. In future work, we will therefore generalize 
this methodology to a wide range of S/N, spatial reso- 
lutions, and spatial sampling in order to investigate the 
effectiveness of our criteria for a more complete range of 



redshifts and observing conditions. 

6. CONCLUSIONS 

We present a simple set of kinematic criteria that can 
be used to distinguish mergers and disks in the SINS sur- 
vey and in similar observations (i.e. with Keck/OSIRIS). 
The reliance of our criteria on dynamical information, 
rather than on surface brightness distributions, takes full 
advantage of the wealth of information provided with in- 
tegral field data. 

We show, via a large set of template galaxies, that 
our criteria can reliably distinguish the majority of ma- 
jor mergers (80%) from disks. When applied to the 
SINS galaxies, this tool provides quantitative support 
for the rotating disk / smooth accretion scenario that 
the interpretation of recent results has suggested (Forster 
Schreiber et al. 2006; Genzel et al. 2006; Birnboim et al. 
2007). In the subset of the SINS systems studied here, we 
quantitatively show that > 50% likely have not had a ma- 
jor merger in their recent history to fuel their rapid star 
formation, providing direct evidence for this scenario. 

Looking forward, the differentiating of disks and merg- 
ers on a galaxy-by-galaxy basis that is now possible is 
useful both immediately and in the future with next- 
generation high-z surveys on 30m class telescopes. As 
an increasing number of galaxies in the high-redshift Uni- 
verse are probed with spatially resolved kinematics, the 
tool presented here can be used to observationally con- 
strain merger fractions, as well as to understand the ef- 
fect of mergers on star formation rates, nuclear activity, 
and growth of structure within proto-galaxies. 
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APPENDIX 

In this appendix, we present the kinemetric analysis on all template and SINS galaxies. These results are summarized 
in Table 2 and are shown for each system in Figures 8 through 14. 



Kinemetry of High-Redshift Galaxies 17 



TABLE 2 
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a Listed here is the 68% confidence interval, as derived from 1000 Monte Carlo realizations, with the three columns indicating the 
lower bound on the confidence interval, the median value, and the upper bound on the confidence interval, respectively. 
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Fig. 8. — Kinemetry of "redshifted" observed template disks at z ~ 2. For each system, the velocity and dispersion fields are shown, 
followed by results from the kinemetric expansion. Overplotted on the velocity and dispersion fields are the major and minor axes of the 
kinemetry ellipses used in the expansion, centered on the continuum center of the system (see text for definition). Because the analysis 
described in Sections 2.2 and 2.3 does not lend itself to straightforward error propagation, we do not include error bars on these figures. 
The most reliable measurement of the kinemetry errors comes from the Monte Carlo realizations, whose results are summarized in Table 2. 
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FlG. 9. — Kinemetry of "redshifted" template mergers at 
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FlG. 10. — Kincmctry of template toy model disks at z ~ 2. Simulation 1 is of intermediate inclination with a centrally-peaked azimuthally 
symmetric light distribution. Simulation 2 is of intermediate inclination and has a light distribution much more extended on one side of 
the galaxy. Simulation 3 is of intermediate inclination, with the light distribution illuminating only one side of the galaxy. Simulation 4 
is nearly edge-on with a centrally-peaked, azimuthally symmetric light distribution. Simulation 5 is nearly edge-on with a varying light 
distribution. 
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Fig. 11. — Kinemetry of cosmologically-simulated disks at z ~ 2. The halo lettering is taken from Naab et al. (2007), who describe the 
evolution of Halos A, C, and E in detail. The snapshots at specific redshifts were selected such that the halos are accreting smoothly, with 
no major mergers, at this point in their histories. 
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Fig. 12. — Kinemetry of cosmologically-simulated mergers at z ~ 2. The halo lettering is taken from Naab et al. (2007), who describe 
the evolution of Halos A, C, and E in detail. The snapshots at specific redshifts were selected such that the halos are undergoing major 
mergers at this point in their histories. 
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Fig. 13. — Kinemetry of SINS galaxies found to be disk-like. 
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Fig. 14. — Kinemetry of SINS galaxies found to be merger-like. 



